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Abstract — LWAl is a new radio telescope operating in the 
frequency range 10-88 MHz, located in central New Mexico. 
The telescope consists of 258 pairs of dipole-type antennas 
whose outputs are individually digitized and formed into beams. 
Simultaneously, signals from all dipoles can be recorded using 
one of the instrument's "all dipoles" modes, facilitating all-sky 
imaging. Notable features of the instrument include high intrinsic 
sensitivity (?s 6 kjy zenith system equivalent flux density), large 
instantaneous bandwidth (up to 78 MHz), and 4 independently- 
steerable beams utilizing digital "true time delay" beamforming. 
This paper summarizes the design of LWAl and its performance 
as determined in commissioning experiments. We describe the 
method currently in use for array calibration, and report on 
measurements of sensitivity and beamwidth. 

Index Terms — Antenna Array, Beamforming, Radio Astron- 
omy. 



I. Introduction 

LWAl ("Long Wavelength Array Station 1"; Figure [T) 
is a new radio telescope operating in the frequency range 
10-88 MHz, collocated with the Very Large Array (VLA; 
107.63° W, 34.07° N) in central New Mexico. The telescope 
consists of an array of 258 pairs of dipole-type antennas whose 
outputs are individually digitized and formed into beams. The 
principal technical characteristics of LWAl are summarized 
in Table U LWAl is so-named because it is envisioned to 
be the first "station" of a 53-station long-baseline aperture 
synthesis imaging array known as the Long Wavelength Array 
(LWA), described in [1], [2]. Although the future of the 
LWA is uncertain, LWAl was completed in Fall 2011 [3] 
and is currently operating under the U.S. National Science 
Foundation's "University Radio Observatories" program. 

Contemporary radio telescopes which are also capable 
of operating in LWAl's 10-88 MHz frequency range in- 
clude GEETEE (35-70 MHz), located in Gauribidanur, India 
[4]; UTR-2 (5-40 MHz), located in the Ukraine [5]; VLA 
(74 MHz) [6]; and LOFAR (10-80 MHz), another new tele- 
scope located in the Netherlands [7], [8]. LWAl and LO- 
FAR are both digital beamforming arrays consisting of large 
numbers of dipole-type antennas and comparable sensitivity; 
however the entire collecting area of LWAl is contained within 
a single station, whereas LOFAR is an aperture synthesis 
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Fig. 1. LWAl. The white cargo container beyond the station array is the 
electronics shelter Visible in the background is the center of the VLA. 



array consisting of many smaller stations distributed over a 
large region of Northern Europe. See [3] for a more detailed 
comparison of these instruments. 

This paper describes LWAl design and performance as 
determined in commissioning experiments. First, in Section HIl 
we provide a brief primer on radio astronomy below 88 MHz 
in order to provide context for subsequent discussion. Sec- 
tion |lll] summarizes the design of LWAl. Section HVl describes 
a simple method for array calibration which we have found 
to be effective. This method uses single-dipole observations 
of strong discrete astronomical sources correlated with ob- 
servations made using an "outrigger" dipole located tens to 
hundreds of wavelengths away. Section [V] reports the results 
of beamforming experiments, and compares the results to those 
predicted in previous LWAl design and simulation studies [1], 
[9]. Conclusions are presented in Section IVll 

II. Background; Radio Astronomy below 88 MHz 

The science applications of compact array beamforming and 
small-aperture imaging below 88 MHz are summarized in [10] 
and [11]; they include the study of pulsars, Jupiter, the Sun, 
and the Earth's ionosphere; studies of the interstellar medium 
including radio recombination lines and electron density varia- 
tions; cosmology through observations of the redshifted 21 cm 
line of neutral hydrogen; and searches for as-yet undetected 
phenomena including radio emission from extrasolar planets, 
astrophysical explosions from a variety of mechanisms, and 
other sources of time-variable or impulsive emission. 

Key issues for antenna and receiver design for radio as- 
tronomical instrumentation operating in this frequency range 
are described in [12], which we briefly summarize here. 
In this frequency range, natural external noise is dominated 
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Beams 


4, independently-steerable 


Polarizations 


Dual linear 


Tunings 


2 center frequencies per beam, 
independently-selectable 


Tuning Range 


24-87 MHz 

(>4:1 sky-noise dominated), 
IU-05 MHz usable 


Bandwidth 


< 16 MHz X 2 tunings X 4 beams 


Spectral Resolution 


Time-domain "voltage" recording; also 
real-time 1024-channel spectrometer 


Beam FWHM 


< 3.2° X [(74 MHz)/j/]^-^ 
(upper bound independent of Z) 


Beam SEFD 
Beam Sensitivity 


~ 6 kjy at Z = 0; depends on pointing, 
celestial coordinates, & frequency; 
see Figure [721 

8 Jy (So-) for 1 s, 16 MHz, Z = 
(inferred from SEFD) 


All-Dipoles Modes 


"TBN": 70 kHz from every dipole, 
continuously 

"TBW": 78 MHz from every dipole, 
in 61 ms "bursts" every 5 min 



Notes: Z is zenith angle, u is frequency. 1 Jy = 10"^'' W m"^ Hz'^. 
FWHM is full-width at half-maximum. SEFD is system equivalent flux density 
(see text). Additional information available at |http : / / Iwal . inf o| and 
http : / / Iwa . unm . edu. 

TABLE I 

LWAl Technical Characteristics. 



by the very bright Galactic synchrotron background, which 
generates antenna temperatures on the order of 10'^ K to 10'' K, 
increasing with decreasing frequency. Contributions from the 
cosmic microwave background (sa 3 K) and other sources 
are also present, but are typically negligible in comparison. 
Man-made noise (other than deliberate radio signals) is also 
negligible for the rural locations at which these instruments are 
typically deployed. Mechanically-steered dishes prevalent at 
higher frequencies are not desirable due to the very large size 
required to achieve acceptably narrow beamwidth. Instead, 
arrays of low-gain dipole-type antennas are used to facilitate 
electronic or digital beamformingfl Individual dipole-type 
antennas typically have impedance bandwidth which is much 
less than the approximately 9:1 bandwidth implied by a tuning 
range of 10-88 MHz; however this is not a limitation as long as 
receiver noise temperature is sufficiently small that the ratio of 
external to internal noise is large after the antenna impedance 
mismatch. A dipole-type antenna combined with a receiver 
having system temperature less than « 500 K is able to achieve 
the best possible (i.e.. Galactic noise-limited) sensitivity over 
a large portion of the 10-88 MHz frequency range. This is 
demonstrated by example in [13] and later in this paper 

External-noise dominance makes the noise measured at 
the output of antenna-receiver channels signicantly correlated 
[9]. The dominance of this correlation over internal noise 
complicates array calibration because the strongest discrete 
astrophysical sources which would otherwise be suitable as 
"incident plane wave" calibrators are orders of magnitude 
weaker than the non-uniform Galactic noise background. The 
method described in Section |IV] of this paper bypasses this 

'The VLA 74 MHz system uses the 27 25-m dishes of the VLA with dipole 
feeds, which yields very low aperture efficiency. This is tolerated as it allows 
reuse of existing high-frequency infrastracture. 




Fig. 2. Intensity of the 74.03 MHz radio sky as measured using LWAl's 
PASI backend (see Section IllTi l, 20;00 local sidereal time (LST). The center 
of this display is the zenith, the perimeter is the hoiizon. North is up. East 
is left. The color scale ranges from ~ 2000 K (dark blue) to ~ 250, 000 K 
(red). The Galactic center ("GC" in this display) is prominent near the bottom 
of the figure, and the Galactic plane extends high into the sky, as indicated by 
the dashed line. Bandwidth; 75 kHz, Integration time: 5 s. "Jup" is Jupiter, 
which is located just below the horizon, and the unlabeled source in the upper 
light is inteiference. The lighter blue regions along the lower left and upper 
right are image noise associated with the point spread function, which has 
not been deconvolved from these images. 

problem by correlating antennas in the array with a distant 
outrigger antenna, which has the effect of suppressing the 
contribution of bright features which are distributed over 
large angular extent - in particular, the Galactic background 
emission which is concentrated along the Galactic plane and 
which is brightest at the Galactic center (see Figures |2] and |3]l. 

Separately from difficulty in calibration, external noise 
correlation significantly desensitizes the beams formed by 
the array [9]. Furthermore, the extent to which a beam is 
desensitized is a function of sidereal time, since the Galactic 
noise intensity varies both spatially and diumally, as is shown 
in Figures |2] and |3] This is particularly frustrating as absolute 
calibration of in-beam flux density is often desired, but now 
depends on both pointing in zenith angle (Z) due to the beam 
pattern variation as a function of Z, and with sidereal time due 
to the spatially- and diurnally-varying Galactic noise intensity 
arriving through sidelobes. 

III. Design 

Antennas. LWAl antennas are grouped into "stands", each 
consisting of a linear-orthogonal pair of antennas, feedpoint- 
mounted electronics, a mast, and a ground screen as shown 
in Figure |4] Each antenna is a wire-grid "bowtie" about 3 m 
long, with arms bent downward at 45° from the feedpoint 
in order to improve pattern uniformity over the sky. The 
feedpoint is located 1.5 m above ground. The ground screen 
is a 3 m X 3 m wire grid with spacing 10 cm x 10 cm 
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Fig. 3. Same as Figure |2] but 14 liours later (10:00 LST). At this time, 
the Galactic center is below the hoiizon and the Galactic plane is low on the 
horizon. 



and wire radius of about 1 mm. The primary purpose of 
the ground screen is to isolate the antenna from the earth 
ground, whose characteristics vary significantly as function 
of moisture content. Zenith values of the effective aperture 
of a single dipole within a stand, multiplied by the loss due 
to mismatch with the lOOfi input impedance of the front 
end electronics and accounting for radiation efficiency due to 
ground loss, are estimated to be 0.25 m^, 8.72 m^, and 2.48 m^ 
for 20 MHz, 38 MHz, and 74 MHz, respectively [9]. Dividing 
out the impedance mismatch loss determined separately [14], 
the associated effective apertures are estimated to be 12.5 m^, 
19.0 m^, and 5.4 m^, respectively. However values for the 
actual in situ antennas are found to vary on the order of 25% 
from antenna to antenna due to mutual coupling. Additional 
details of the design and analysis of these antennas can be 
found in [9] and [14]. 

Front end electronics. Each dipole is terminated into a pair 
of commercial InGaP HBT MMIC amplifiers (Mini-Circuits 
GALI-74) in a differential configuration, presenting a lOOil 
balanced load to the antenna. This is followed by a passive 
balun which produces a 5051 single-ended signal suitable for 
transmission over coaxial cable, plus some additional gain 
to overcome cable loss. The total gain, noise temperature, 
and input 1 dB compression point of the resulting front 
end electronics units are approximately 36 dB, 300 K, and 
— 18 dBm respectively, and are approximately independent 
of frequency over 10-88 MHz [14]. The gain and noise 
temperature of the feedpoint electronics are such that they 
dominate the noise temperature of the complete receiver chain, 
which is much less than the antenna temperature as is desirable 
(See Section HH). The 1 dB compression point has been found 
to be satisfactory at the LWAl site. Although higher 1 dB 
compression would be better, this would be difficult to achieve 




Fig. 4. LWAl antenna stands. Front end electronics are enclosed in the white 
box at the feedpoint. Signals exit through coaxial cables inside the mast. Each 
stand is about 1.5 m high. 



without compromising noise temperature. Additional infor- 
mation on design requirements for the front end electronics 
is available in [15]. The curve labeled "Sky" in Figure |5] 
shows the noise temperature measured by an LWAl receiver 
in normal operation (i.e., dipole attached), after calibration 
to account for the known gain of the electronics following 
the antenna terminals. Thus, this is an estimate of the system 
temperature. Also shown is the same measurement with the 
dipole terminals short-circuited, which is assumed to zero 
the noise delivered to the front end electronics unit without 
significantly changing its behavior, which is consistent with 
the findings of laboratory experiments. These result confirm 
an internal noise temperature of about 300 K, 4:1 external 
noise dominance over 24-87 MHzIl and negligable level of 
intermodulation. 

Array Geometry. LWAl consists of 256 antenna stands 
(512 antennas) within a 100 m (East- West) x 110 m (North- 
South) elliptical footprint, plus two stands (4 antennas) which 
lie outside this footprint (the "outriggers"). The arrangement 
of stands is shown in Figure |6] The station diameter and 
number of stands per station were originally determined from 
an analysis of requirements for the LWA aperture synthesis 
imaging array, as detailed in [1], [2]. However, these choices 
are also appropriate for the present single-station instrument, 
as is demonstrated in this paper. This choice of station aperture 
and number of stands results in a mean spacing between 
stands of about 5.4 m, which is G.36A and 1.44A at 20 MHz 
and 80 MHz respectively. To suppress aliasing, antennas are 
arranged in a pseudo-random fashion, with a minimum spacing 
constraint of 5 m in order to facilitate maintenance. The 
elongation of the station aperture in the North-South direction 
improves main lobe symmetry for pointing towards lower 
declinations, including the Galactic Center, which transits at 
Z « 63° as seen from the site. 

-It should be noted that the Galactic background-dominated antenna tem- 
perature varies diurnally over a range of about 35% due to the rotation of the 
Earth (see [13] for an illustration); as a result the external noise dominance 
varies slightly over a 24-hour period. 
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Fig. 5. Power spectral density measured by an LWAl receiver, calibrated 
to the antenna terminals. Top curve: Result expressed as system temperature 
(spikes are human-generated signals); Bottom cun'e: Same measurement made 
with a short circuit termination at the input, which provides an estimate of the 
internal contribution to the system temperature. Spectral resolution: 6 kHz. 
Integration time: 10 s. Early afternoon local time. 
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Fig. 6. Arrangement of stands in the LWAl an'ay. The minimum distance 
between any two masts is 5 m (0.33A, 0.63A, and 1.23A at 20 MHz, 38 MHz, 
and 74 MHz, respectively). All dipoles are aligned North-South and East- West. 
Outiigger stands are not shown. 



Cable System. Connections between the front end electron- 
ics and the electronics shelter are by coaxial cables having 
lengths between 43 m and 149 m (excluding outriggers). 
These cables have loss of about —0.11 dB/m at 88 MHz, 
and dispersive delays (that is, frequency-dependent delays in 
addition to the delay implied by the velocity factor) given by 
[16], [17]: 



^100 my VlO MHz/ 

where I is length and ly is frequency. Thus the signals arriving 
at the electronics shelter experience unequal delays, losses, 
and dispersion. These can be corrected either "in line" in 




Fig. 7. View from within the LWAl array looking toward the equipment 
shelter 



the station's digital processor (see below) or, for the all- 
dipoles modes, as a post-processing step. In the current real- 
time beamforming implementation, the non-uniform dispersive 
delays are compensated for the center frequency of the highest- 
frequency tuning in a beam; thus there is some error over the 
bandwidth of the beam, and additional dispersion error for the 
lower-frequency tuning in the same beam. 

Receivers. Upon arrival in the shelter (Figure |7), the sig- 
nal from every antenna is processed by a direct-sampling 
receiver comprised of an analog section consisting of only 
gain and filtering, a 12-bit analog-to-digital converter (A/D) 
which samples 196 million samples per second (MSPS), and 
subsequent digital processing to form beams and tune within 
the digital passband. Digitization using fewer than 12 bits 
would be sufficient [15], but the present design eliminates 
the need to implement gain control in the analog receivers 
and provides generous headroom to accommodate interference 
when it becomes anomalously large. The choice of 196 MSPS 
ensures that strong radio frequency interference (RFI) from 
the 88-108 MHz FM broadcast band (see Figure |5]l aliases 
onto itself, with no possibility of obscuring spectrum below 
88 MHz. To accommodate the various uncertainties in the RFI 
environment, analog receivers can be electronically reconfig- 
ured between three modes: A full-bandwidth (10-88 MHz) 
uniform-gain mode; a full-bandwidth dual-gain mode in which 
frequencies below about 35 MHz can be attenuated using a 
shelf filteiH and a 28-54 MHz mode, which serves as a last 
line of defense should RFI above and/or below this range 
become persistently linearity-limiting. In addition, the total 
gain in each mode can be adjusted over a 60 dB range in 2 dB 
steps, allowing fine adjustments to optimize the sensitivity- 
linearity tradeoff. Use of the 28-54 MHz mode has not been 
required to date. 

Digital Beamforming. A detailed description of the LWAl 
digital processor is provided in [18] and is summarized here. 
Beams are formed using a time-domain delay-and-sum archi- 
tecture. Delays are implemented in two stages: An integer- 

^A "shelf filter" is a filter which has two adjacent passbands, with one 
passband (the shelf) having higher attenuation than the other 
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sample "coarse" delay is applied using a first-in first-out 
(FIFO) buffer operating on the A/D output samples, followed 
by a 28-tap finite impulse response (FIR) filter that implements 
an all-pass subsample delay. The filter coefficients can also 
be specified by the user, allowing the implementation of 
beams with custom shapes and nulls. The delay-processed 
signals are added to the signals from other antennas processed 
similarly to form beams. Four dual-polarization beams are 
constructed in this fashion, each fully-independent and capable 
of pointing anywhere in the sky. Each beam is subsequently 
converted to two independent "tunings" of up to 16 MHz 
bandwidth (4-bits "I" + 4-bits "Q" up to 19.6 MSPS) each, 
with each tuning having a center frequency independently- 
selectable from the range 10-88 MHz. Both polarizations 
and both tunings of a beam emerge as a single stream of 
user datagram protocol (UDP) packets on 10 Gb/s ethernet. 
Thus there are four ethernet output cables, with each one 
representing two center frequencies from a particular pointing 
on the sky. The maximum data rate (ignoring protocol bits) on 
each ethernet cable carrying beam data is therefore 19.6 MSPS 
X 8 bits/sample x 2 polarizations x 2 tunings = 627.2 Mb/s. 

All-Sky Modes. Simultaneously with beamforming, LWAl 
is able to coherently capture and record the output of all its 
A/Ds, where each A/D corresponds to one antenna. This can be 
done in two distinct modes. The "transient buffer - wideband" 
(TBW) mode allows the raw (« 78 MHz) 12-bit output of 
the A/Ds to be collected in bursts of 61 ms at a time, and 
« 5 minutes is required to write out the captured samples. 
The "transient buffer - narrowband" (TEN) mode, in contrast, 
allows a single tuning of w 70 kHz bandwidth to be recorded 
continuously for up to 20 hours. These two modes share the 
same 10 Gb/s ethernet output from the digital processor, and 
thus are mutually exclusive. However, the TBW/TBN output 
is distinct from the four beam outputs and runs simultaneously 
with all four beams. 

Data Recorders/Spectrometers. The Umited data rate of the 
internet connection from the LWAl site makes data transfer 
from the site impractical for observations longer than a few 
minutes. Instead, each beam output and the TBW/TBN output 
is connected to a dedicated data recorder (DR). A DR is a 
computer that records the UDP packets to a "DR storage 
unit" (DRSU). Currently, a DRSU consists of five 2 TB drives 
(10 TB total) in a lU rack-mountable chassis, configured as 
an eSATA disk array. Each DR can host 2 DRSUs. At the 
maximum beam bandwidth, each DRSU has a capacity of 
about 30 hours of observation. The data saved is fully coherent 
and Nyquist-sampled. Alternatively, the DRs can be used in 
"spectrometer mode" in which they continuously compute 32- 
channeQ fast Fourier transforms (FFTs) on the incoming beam 
data independently for each beam and tuning, and then time- 
average the EFT output bins to the desired time resolution. 
This results in a dramatic reduction in data volume, but is not 
suitable for all observing projects. 

PASI. The "Prototype All-Sky Imager" (PASI) is a software- 
defined correlator/imager currently operational at LWAl using 
the TBN data stream. It consists of a cluster of 4 server- 

''increased to 1024 channels since this paper was originally submitted. 



class computers with Nehalem multicore processors inter- 
connected by an Infiniband switch. PASI images nearly the 
whole sky in the Stokes / and V parameters many times 
per minute, continuously and in real time, with an average 
duty cycle of better than 95%. (The V parameter is useful 
for discriminating circularly-polarized emission such as that 
from Jupiter) PASI does this by cross-correlating the dipole 
data streams, producing a sampling of "visibilities" within the 
station aperture. These visibilities are then transformed into 
sky images using the NRAO's Common Astronomy Software 
Applications (CAS A) software libraryjf] Figures |2] and [3] were 
obtained from PASI. 

IV. Array Calibration 

In this section we describe the technique currently in use 
to calibrate the array for real-time beamforming. Since LWAl 
uses delay-and-sum beamforming, the problem is to determine 
the set of delays which, when applied to the dipole (A/D) 
outputs, results in a beam with maximum directivity in the 
desired direction, subject to no other constraints. (In [9], this 
is referred to as "simple" beamforming. Other approaches are 
possible but are not considered in this paper) In principle these 
delays can be estimated a priori, since the relevant design 
values (in particular, cable lengths) are known in advance. 
In practice this does not work well due to errors in presumed 
cable lengths and the accumulation of smaller errors associated 
with distribution of signals within equipment racks. Thus, 
delays must be measured using external stimulus signal(s) 
while the instrument is in operation. 

The approach used here is to decompose the problem into a 
set of narrowband calibration problems which are solved using 
data collected using LWAl's TBN mode, and then to extract 
delays by fitting a model of the presumed cable response 
(including dispersion) to the measured phase responses. The 
narrowband decomposition is justified by the fact that the 
maximum time-of-flight between any two antennas in the 
station array is < 367 ns, which is much less than the TBN 
inverse bandwidth (70 kHz)~^ w 14 /is. As described below, 
the outriggers play an important role in the calibration; the 
maximum time-of-flight between any antenna in the station 
array and the outrigger (maximum separation « 390 m) is 
< 1.3 fis. 

The narrowband procedure is described in rigorous mathe- 
matical detail in [19], and summarized here. Each narrowband 
calibration relies on the ability to identify the response due 
to a point source in the array output. This is confounded by 
the problem that there are typically multiple bright sources 
present in the sky, and (as explained in Section|Il]i the problem 
that Galactic noise dominates the system temperature and thus 
appears in the data as a bright, non-uniformly distributed 
source. To suppress the effect of distributed features (including 
the Galactic center and Galactic plane) as well as discrete 
sources other than the source of interest, we use a "fringe rate 
filtering" technique. Fringe rate filtering is essentially a simple 
narrowband version of the delay / delay-rate filtering technique 
of Parsons & Backer (2009) [20]. In this technique, each 

^http://casa.nrao.edu/ 
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Fig. 9. Fringe rate spectram corresponding to Figure |8] This result is 
dominated by two bright sources (these happen to be the radio galaxy 
Cygnus A (Cyg A) and the supernova remnant Cassiopeia A (Cas A)) and 
a "DC" term. The DC term is the "all sky" contribution predicted in [9], 
and stands as additional evidence that the obsereation is strongly external 
noise-dominated. Many other sources ai'e present, but ai'e orders of magnitude 
weaker The frequency resolution is 279 /iHz. 



Fig. 8. Fringes between a dipole on the far West side of the array and an 
outrigger dipole located 389 m to the East. Top: Magnitude. Bottom: Phase. 
Each point represents 10 s integration. Start time is 18:50 LST. 



dipole in the array is correlated with a corresponding outrigger 
dipole. The electrically-large spacing makes the correlation 
relatively insensitive to spatial structure in the sky noise 
intensity having large angular scales (see [9] for examples.) 
The contributions to the correlation due to individual discrete 
sources exhibit time-varying phase due to the apparent rotation 
of the sky; in astronomical parlance, these source-dominated 
correlations are referred to as "fringes"; see Figure 0for an 
example. The rate of phase rotation in the fringes depends 
on the positions of the antennas being correlated relative 
to the direction to the sources. A time-to-frequency Fourier 
transform of the correlation over an interval greater than the 
reciprocal of the smallest source-specific fringe rate yields a 
"fringe rate spectrum", in which discrete sources are apparent 
as localized components; an example is shown in Figure |9] 
"Fringe rate filtering" refers to the process of selecting just 
one of these components, suppressing all others (which can 
be done by matched filtering ("fringe stopping") in the time 
domain or by excision in the fringe rate domain); the resulting 
correlation then represents the response of just one source. 
This of course presumes that sources are sufficiently separated 
in fringe rate: This is one reason for the use of outriggers, 
since fringe rate is proportional to antenna separation. Further 
improvement is possible by proper scheduling of the associated 
observations; that is, choosing times in which suitable sources 
have sufficiently different fringe rates. For completeness, the 
positions and strengths of the sources used in the example 
shown in Figures |8] and |9] are given in Table |II] 

This procedure is repeated for every antenna in the ar- 
ray, yielding a set of complex-valued coefficients associated 
with one direction at one frequency. One then iterates over 



Source 


RA 


Dec 


Az 


Z 


Flux density 


Cyg A 

Cas A 


19h 59' 
23h 23' 


-1-40° 44' 
-|-58°48' 


46°08' 
38°29' 


10° 11' 
46° 27' 


17.06 kJy 
17.13 kJy 



Right a.scension (RA) and declination (Dec) are given for tlie J2000 epoeli. 
Azimutli (Az) and zenith angle {Z] are indicated for the midpoint of the observation. 

The flux densities are the 73.8 MHz values given in LWA Memo 155 [21], scaled to 74.03 MHz using spectral indices of 
— 0.58 and -0.72 for Cyg A and Cas A respectively, and Cas A is further scaled by -O.TS^/yr [22], 



TABLE 11 

Astronomical source data pertaining to Figures[8]and[9] 



frequency. There is currently no iteration over direction; 
that is, only one direction is considered. Different directions 
give different coefficients, for two reasons: First, because the 
geometrical delays are different; however these are easily 
equalized. The second reason is because the antenna patterns 
are made unequal by the effects of mutual coupling. The 
effects of mutual coupling on the antenna patterns is difficult 
to know precisely. However since our goal is to compute 
(broadband) delays, as opposed to (narrowband) phases and 
magnitudes, calibration in a single direction is sufficient, as 
will be demonstrated in the next section. The entire procedure 
requires on the order of a week, typically observing over the 
same one-hour LST range at a different frequency each day. 

V. Beamforming Performance 

A. Multi-frequency Drift Scan at Z — 6.7° 

Figures [Tol and [TTI show a single transit drift scan of Cyg A; 
i.e., the output of a fixed beam pointed at the location in 
the sky at which Cyg A achieves upper culmination: Z = 
6.7°, 0° north azimuth. Cyg A is useful as a test source 
both because it is the strongest time-invariant and unresolved 
astronomical source visible from the LWAl site, and also 
because it transits close to the zenith. The origin in both 
plots corresponds to the time at which Cyg A is expected 
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0.4 - 



-60 -40 -20 20 40 

Time [min] 

Fig. 10. Simultaneous drift scans of Cyg A at 85.00, 74.03, 62.90, 52.00, 
and 37.90 MHz. Scans can be identified by peak width, which increases 
with decreasing frequency. All scans are normalized to a common maximum. 
Bandwidth: 211 kHz, Integration time: 786 ms. Single (N-S) polarization, no 
RFI mitigation. 




0.4 - 



-60 -40 -20 20 40 

Time [min] 

Fig. 11. Simultaneous drift scans of Cyg A at 28.80 and 20.50 MHz. Scan 
magnitudes are normalized. Bandwidth: 211 kHz, Integration time: 786 ms, 
Single (N-S) polarization, no RFI mitigation. 



to peak; Figure |2] shows the sky at this time. The results for 
a narrow bandwidth at 8 distinct center frequencies ranging 
from 20.5 through 85.6 MHz are shown. The large peak in 
each curve is Cyg A; the small peak to the right in some curves 
is the Galactic plane passing through the beam a short time 
later (apparent only at higher frequencies and unresolved in 
lower-frequency curves). Note that Cyg A is clearly detectable 
at all frequencies, demonstrating the very wide instantaneous 
bandwidth of the instrument. Also apparent is the broadening 
of the beamwidth with decreasing frequency, as expected. The 
lowest two frequencies (28.80 and 20.50 MHz) are shown in 
a separate figure for clarity; both frequencies are apparently 
strongly affected by ionospheric scintillation, and 20.05 MHz 
is strongly affected by RFI. 

The drift scans for each frequency in Figure [TO] and [TT] are 
normalized to the same peak value to facilitate comparison 
of beam width. Using the method described in Appendix |A] 



we find full-width at half-maximum (FWHM) beamwidths of 
2.8°, 3.5°, 4.6°, 5.8°, 7.5°, 8.0° and 10.9° for 85.00, 74.03, 
62.90, 52.00, 37.90, 28.80, and 20.50 MHz respectively. Note 
that these measurements are in planes of constant declination; 
i.e., not planes of constant azimuth nor constant elevation. The 
expected zenith FWHM of a uniformly-illuminated circular 
aperture having the same minimum dimensions as LWAl is 
2.3° at 74 MHz; thus the measurements are roughly 20-60% 
larger than predicted by the above model. This is due in part to 
the proximity of Cyg A to the Galactic plane, and the resulting 
"confusion" of the Cyg A emission with diffuse Galactic 
emission filling the beam. This is evident by comparison with 
the same measurement performed using the source 3C123, 
which lies further from the Galactic plane but transits at a 
similar zenith angle (Z 4.4°). For 3C123 we obtain FWHM 
= 2.7° at 74.03 MHz, which is quite close to the expected 
value. Due to this limitation in the measurement technique, 
the FWHM values reported here should be considered upper 
bounds as opposed to estimates of the actual values. 

We now consider beam sensitivity. Because LWAl is 
strongly Galactic noise-limited, beam directivity is not a 
reUable metric of sensitivity. Instead, we consider sensitivity 
in terms of system equivalent flux density (SEFD), which 
is defined as the strength of an unresolved (point) source 
that doubles the power at the output of the beam relative 
to the value in the absence of the source. Because the flux 
density of Cyg A is known (Table beam SEFD can be 
determined directly from drift scans using the method de- 
scribed in Appendix IbI This method yields SEFD = 16.1 kJy 
for 74.03 MHz. To obtain the source flux S at some other 
frequency v, we apply the known spectral index —0.58 for 
Cyg A [23]: 

/ ^ -0.58 

S{y)^S{v^)i^] , (2) 

where i/q is a frequency at which the flux is already known. 
From this we estimate SEFD = 11.0, 21.5, 28.8, 21.5, 18.2, and 
47.0 kJy for 85.00, 62.90, 52.00, 37.90, 28.80, and 20.50 MHz 
respectively. Due to the variables described in Section |II] it 
is difficult to know even within lO's of percent what values 
to expect; however it is shown in the next section that these 
results are consistent with the results of simulations in previous 
work. It should also be emphasized that the flux density of 
radio sources (including Cyg A and sources discussed in the 
next section) is known only approximately below 74 MHz, 
and that use of the 74 MHz flux density and spectral index to 
calculate flux density at lower frequencies is another source 
of error. 

B. Sensitivity & Beamwidth vs. Frequency & Elevation 

The same procedures described in the previous section have 
been used to determine sensitivity at lower elevations using 
Cyg A and the strong sources Tau A, Vir A, 3C123, and 
3C348. The results are summarized in Figure [T2] Also shown 
in this figure are predictions originally shown in Fig. 8 of 
[9]. Note that SEFD depends both on pointing with respect to 
celestial coordinates as well as sidereal time (for the reasons 
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Fig. 12. Sensitivity (SEFD) vs. Z obtained from diift scans. The numbers 
used as data markers indicate the frequency rounded to the nearest MHz. 
Markers in red italic font represent transit drift scans of the sources 3C123 
(Z = 4.4°), Tau A (12.1°), Vir A (21.8°), and 3C348 (29.1°); all others 
are Cyg A. The curves are predictions from Fig. 8 of [9] for (bottom to top) 
74, 38, and 20 MHz. The use of both polarizations is assumed. 

explained in Section|II]l; whereas the predictions of [9] assume 
a uniformly bright sky. Thus, precise agreement cannot be 
expected. Nevertheless, the levels and trends with Z and 
frequency appear to be consistent with the predictions. We 
do note however that the agreement for 38 MHz appears to 
be somewhat worse than the agreement at other frequencies. 
We do not currently have an explanation for this. 

Figure [T3] shows a summary of FWHM measurements from 
drift scans. Beamwidth is difficult to measure at large Z 
because the drift scan peaks become simultaneously broad 
and weak; thus the data shown is limited \a Z < 45°. Note 
that Cyg A measurements generally indicate FWHM higher 
than expected, whereas measurements using other sources are 
closer to expectations; again this is due in some part to the 
inability of the method of Appendix |A] to account for the 
excess emission associated with the nearby Galactic plane. 
Also shown in Figure [13] is an empirical upper bound based 
on all measurements (regardless of Z) considered to date: 

FWHM < 3.2° X (3) 

We intend to improve on this bound in future measurements 
by using interferometric methods, as discussed in Section |VI] 



VI. Conclusions 

This paper described the design of the LWAl radio tele- 
scope, array calibration technique, and the results of com- 
missioning experiments to confirm beamforming performance. 
The "single point" delay calibration technique currently in use 
at LWAl was described in Section |IV] In this technique, the 
delays are estimated from phase estimates over a range of fre- 
quencies using narrowband "all-dipoles" (TBN) observations 
of the sky. This technique uses a fringe rate filtering technique 
which is robust to the presence of multiple bright sources 



Fig. 13. Beamwidth (FWHM) vs. frequency obtained from drift scans. 
Markers indicate Z: "o", 6.7°; "+", 25°; squai'e, 45°; all using Cyg A. 
"x" are 74.03 MHz results using transit drift scans from the same additional 
sources identified in the caption of Figure 1121 all at 74.03 MHz. The solid 
line is the upper bound described in the text. 

including both unresolved and spatially-extended sources. Real 
time delay-and-sum beamforming was demonstrated in Sec- 
tion [V] Sensitivity findings are summarized in Figure [12] and 
were shown to be consistent with predictions made in [9]. 
Beamwidth findings are summarized in Figure [13] Due to the 
limitations of the drift scan measurement technique employed 
in this paper, we can currently only upper-bound our estimates 
of beamwidth. In future measurements we intend to cross- 
correlate drift scan beams with outrigger dipoles, which is 
expected to suppress the contribution of diffuse Galactic emis- 
sion and result in more accurate measurement of beamwidth. 
The same technique will enable characterization of sidelobes, 
which is not possible with the drift scan technique. We also 
plan to conduct measurements to characterize polarization 
performance. The results presented here combined with early 
science results described in [3] confirm that LWAl perfor- 
mance is consistent with expectations and that the instrument 
is ready for science operations. 

Appendix A 
Estimation of Beamwidth from a Drift Scan 

The peak-to-baseline ratios achieved in drift scans such as 
those shown in Figures [10] and [TT] are too low to facilitate 
calculation of beam FWHM by direct inspection. However it 
is possible to estimate FWHM by modeling the drift scan as 
the sum of a Gaussian function representing the power pattern 
of the main lobe plus a constant noise floor, solving for the 
coefficient of time in the exponent of the Gaussian function, 
and then calculating the associated width of the Gaussian 
function alone. In this appendix we first derive the result, and 
then we justify the use of the Gaussian approximation. 

The drift scan is modeled as: 

y{t) = A + Be""^*' (4) 

where A is the magnitude of the constant noise baseline, B 
is the peak magnitude of the Gaussian function modeling the 
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main lobe, and t is time. Evaluating this expression at the time 
t = 0, taken to be the time at which the peak occurs, and time 
ti, a short time later, we have: 



Solving for 7: 



7 



yiO) = A + B 
y{ti) = A + Se"''*! . 

ln(y(0)^A)-lnfa(ti)-A) 



(5) 
(6) 

(7) 



FWHM corresponds to the time interval between the half- 
maximum points of the Gaussian function, which is easily 
found to be 



FWHM 




(8) 



We now show that the Gaussian model is a reasonable 
approximation to the exact result for the zenith-pointed beam 
due to a uniformly-weighted circular aperture. Note the LWAl 
aperture is not exactly circular, nor is it exactly uniformly- 
weighted (due to variations in antenna responses due to mutual 
coupling); our purpose here is simply to demonstrate that 
the Gaussian model reasonably represents the kind of beam 
that we expect to see from this type of antenna system. The 
normalized (i.e., maximum magnitude = 1) pattern function 
of a uniformly-weighted circular aperture lying in the plane 
of the ground is [24]: 



F{Z) = 2 



J I {j3a sill Z) 
Pa sin Z 



(9) 



where a is the radius of the aperture, /3 = 27r/A, and A is 
wavelength. Thus the normalized power pattern is F^{Z). Let 
us assume that the main lobe of F'^{Z) can be modeled as 
a Gaussian function y{Z/c), where c is the angular rate of 
drift. For purposes of demonstration let us assume a/A = 12.3 
(100 m diameter at 74 MHz). We take A = 0, y{Q) = 1, 
y[t^ = Zi/c) = F'^{Zi) = 0.9, and then use Equation |9] to 
find cti = 0.475°. Figure fT4l compares F^{Z) to the derived 
model y{Z/c) with 7 obtained from Equation|7] Note that the 
model error at the half-maximum point is very small, justifying 
the use of the Gaussian model for this purpose. 

Finally, we note that the model has not been demonstrated 
to be appropriate for beam pointings far from the zenith. 
However, we note that we do not have a "true" model available 
for this case either, since the patterns of individual antennas 
combined with mutual coupling make this difficult to know. 
This caveat should be considered when evaluating results 
computed for low-elevation pointings. 

Appendix B 
Calculation of SEED from Drift Scans 

When the flux density of the source is known, beam system 
equivalent flux density (SEED) can be estimated from drift 
scans such as those shown in Figures [TO] and [TT] A derivation 
follows. The power density P, captured in one polarization of 




1 1.5 

Zenith Angle [deg] 



Fig. 14. Comparison of the Gaussian main lobe model y{cZ) (blue/dashed) 
to the main lobe power pattern F'^{Z) (red/solid) for a/A = 12.3 . 



the beam associated with an unpolarized source having flux 
density S is given by 



Fc — S Ap 
2 



(10) 



where A^ is the effective aperture. Similarly, the power density 
N captured in one polarization of the beam associated with 
the total available noise is given by 



(11) 



where k is Boltzmann's constant and Tgys is the system noise 
temperature, assuming this to be external noise-dominated. 
Thus the signal-to-noise ratio at the beam output is 



N 



SA, 



(12) 



sys 



In these terms, the SEED for this polarization is defined as the 
value of S which results in Pg/N = 1; thus: 



SEED 



(13) 



Let Pi = Pg + N he the power density measured at the peak 
of the drift scan, and let Pq = N he the power density in 
the absence of the source. Note Pq can be estimated from the 
approximately-constant noise baseline on either side of the 
peak. Using the above definitions, we find 



Solving for SEED: 



2 sys 



s 



2'^-^ sys 



SEED 



1 



SEED = S 



Pn 



1 



(14) 



(15) 



Thus, given the source flux density, the SEED can be obtained 
from Pi and Pq, which themselves can be read from a single- 
polarization drift scan. 
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